matlab boundary value problem